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2 Summary 

The Riemann C-function is cr > 1, or 

its analytic continuation for ct < 1. It was originally 
studied by Euler on the integers and extended to the 
complex plane by Riemann for the purpose of inves¬ 
tigating the properties of prime numbers. Today, it 
finds application in statistical physics and quantum 
field theory [2]. 

An inherent symmetry underlies the dependence 
of the partial sums on n. This symmetry 

relates individual stepji] n < Up = to 

collections of adjacent steps n > Up. It defines two 
particular points in the complex plane: 

A center point conjugate to itself in a sense 
described below, with an axis of bilateral 
symmetry passing through it at a specific 
angle. 

A point conjugate to the origin, beyond 
which the sum of steps either converges or 
approaches infinity, depending on the value 
of cr. 

^ [x] denotes the integer floor of x 


Forward step angle differences S0n = —t/{n + 
1/2) -I- ..., = tl{n + 1/2)^ -f ..., etc. are de¬ 

fined between adjacent steps. At the point referred 
to here as Up, defined by the partial sum 
the second angle difference equals 27r, equiv¬ 

alent to zero. As n increase beyond np, the second 
angle variation S^9n decreases monotonically from 27r 
to 0, while first angle variation S9n decreases through 
Up regions bounded by steps whose first angle differ¬ 
ences are odd multiples of tt. Each of the Up regions 
[rip/(n-l-1)] < N < [np/27r(n— 1)] is conjugate to the 
single initial step n. There are approximately n^/n^ 
steps in the region conjugate to n. Because of trun¬ 
cation to step numbers n G P, the precise boundaries 
of conjugate regions are not exactly coincident with 
centers of steps. 

When Up » 1, angles and lengths of the steps im¬ 
mediately adjacent to rip are nearly constant, form¬ 
ing a polygon tangent to step Up. P{s), the center 
of symmetry, is recognized here to be the center of 
this feature, laterally displaced from Up. The spe¬ 
cific form of the conjugacy is that each initial step 
n“®, transformed to and multiplied by a com¬ 
plex factor independent of n, equals the sum of steps 
in its conjugate region. The polar form of the com¬ 
plex factor is Q{s) = (t/27r)^/^“'^e^*®, where 0(s) is 
approximately the sum of the angles of steps Up and 
rip+i. To the next order of approximation, end cor¬ 
rections A„ for the conjugate regions give the precise 
points conjugate to the ends of step n, assumed to 
correspond to the geometric center of the constella¬ 
tion of steps having 59 ~ the nth odd multiple of tt at 
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Nn- The lengths of a conjugate region and its initial 
step n are equal when \Q\ = 1, which only occurs for 
cr = 1/2. Q{s) 7 ^ 0 on the critical strip. The sum of 
the angles of each initial step and its summed conju¬ 
gate region is 20(s) and the axis of symmetry is at the 
angle 0 ± 7r/2. Angle differences of adjacent regions 
conjugate to n, n -|- 1 can be summed immediately, 
and are seen to equal S9n- The sum of steps in a con¬ 
jugate region requires explicit summation, which can 
be approximated by a Jacobi 0-function for a large 
number of conjugate steps. This is the case for t 1 
and n <C Up. 

The description of C(s) in terms of symmetry is 
shown in subsequent sections to agree with analytic 
results. The continuation of C(s) for a < 1 corre¬ 
sponds to the point 0'{s) conjugate to the origin, 
using Riemann’s general definition. As a specific ex¬ 
ample, O' —>■ C(s) for cr < 1 is consistent with the 
Euler-Maclaurin procedure. The functional equation 
C(s) = (5(s)C(l — s) follows from the symmetry. The 
central symmetry point P{s), which is conjugate to 
itself, is related directly the Riemann-Siegel equation; 
the lowest order geometric construction reproduces 
Riemann’s first order of correction. Other specific 
results obtained by Riemann and his successors are 
addressed, including the number of zeros of /'(ct -I- it) 
for t < T as a function of 0(s) and their relation to 
the critical line a = 1/2. 

In this work, detailed derivations and figures are 
omitted in favor of brevity. A more complete discus¬ 
sion, including bounds, generalized C functions and 
the special significance of prime steps, is found in [7]. 


3 A Review of the Traditional 
Analysis 

Mathematicians of the early 19th century obtained 
general expressions for functions previously defined 
only on the integers, and extended them to the com¬ 
plex planeH Euler derived an expression for the fac¬ 
torial in terms of the definite integral familiar to all 
mathematics students 

^See, for example, [J] Chap. 12]. These topics are chosen 
both to illustrate the history and to introduce specific functions 
encountered below. 


pOO 

r(a:) = / —>■ (ic — 1)! for xSP (1) 

Jo 

which Gauss expressed alternately as 


r(a:) = lim 


1 • 2 • ... (iV - 1) 


N^oo (x + l)(x -b 2) . . . (x -b iV — 1) 


( 2 ) 


Euler’s integral arises naturally in analysis, and 
Gauss’ form is useful to define connections between 
T-functions. For example, the functional equation 
r(x -b 1) = xr(x) is supplemented by sirnrs = 
7rs/r(s)r(l - s) and r(s) = 2r((s -b l)/2)r(s/2) . 

A discrete series can often be represented by the ex¬ 
pansion of a continuous generating function^] this 
facilitates integrations in the complex plane. An el¬ 
ementary application from statistical physics is the 
integral of the Planck distribution for photon^ 


x^dx 
e^ -I 


dx, 


(3) 


Expansion of the integrand gives the series under 
discussion here with argument s = 4: 


poo ^ OO 

/ x3^e""dx = r(4)^n-4^r(4)C(4) (4) 

Jo 1 1 

Riemann generalized this to the complex plane by 
replacing the numerator x^ of this example by {—x)" 
and using a contour from oo above the positive real 
axis, circling the origin in the positive sense at radius 
S, followed by a return to c» below the positive real 
axis. He showed that this procedure yields analytic 
continuation for tr < l[3l Sec. 1.4]. 

Euler had earlier expressed the fundamental 
Dirichlet series as a product of terms involving the 
prime numbers (his “Golden Rule”) 



®The numerator x^dx arises from the dimensionless form of 
the photon energy hu times {hujc)^ d{hi/jc) for a volume ele¬ 
ment of momentum space, and the denominator represents the 
quantum distribution function for bosons, where x = hv/kT. 
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Riemann used the poles resulting from the loga¬ 
rithm of this product as the starting point of his anal¬ 
ysis leading to a relation between zeros of C and the 
primes, as well as an expression for n(a:), the number 
of primes less than or equal to a; [31 Sec. 1.18] 

n(a;) = Li{x) — ^ Li{x^) -f ... (6) 

zerosp 

His first proof of the functional equation of the 
C-function under the s —>■ 1 — s transformation [31 
Sec. 1.6] 

C(s) = ^^C(l-s) (7) 

r(s) 

used the contour integration above, in addition to dis¬ 
joint circles about the poles of the generating function 
on the imaginary axis in the negative sense. The two 
contributions add to zero, giving Equation [T] It will 
be seen that this is a detailed relation, applying not 
only globally to the sum over all steps and poles, but 
to the individual pairs of terms as well. 

Jacobi’s 0-functions include 

CXD OO 

G{u) = and T(u) = Y e""'™ (8) 

— CO 1 

which Riemann used for the second proof of his func¬ 
tional equation [31 Sec. 1.7]. He then applied Gauss’ 
T-function identities to express a function ^(s) which 
is invariant under s —>■ 1 — s as 

ax) = {s-l)T{l + s/2)7r-^/\{s) (9) 

The history of the use of these integrations and 
inversion transforms to examine the prime numbers is 
described in detail by Edwards [3] . General methods 
of analytic number theory are discussed by Apostol 

a- 

4 Angles and Conjugate Re¬ 
gions 

For extremely large argument t and n <C I/tt, the first 
angle differences contain many multiples of 27r. Al¬ 
though step angles vary monotonically as —tlog{n), 


their modulus with respect to 27r appears to vary ir¬ 
regularly when S9 ^ 27r. For one particular num¬ 
ber Up of initial steps, there are Up subsequent steps 
having 69 equal to an even multiple of tt. Solving 
Up = 69n^l2'n gives Up = [y^t/27r- l/2]0 For rip > 1 
its adjacent steps approximate the “pendant”, a reg¬ 
ular polygon with slowly varying lengths and 69, a 
distinct feature in the step plot. See Figure [21 For 
all steps n > Up, 69n decreases smoothly through odd 
and even multiples of tt, giving Up regions conjugate 
to the Up initial steps. The center step of the region 
conjugate to step n is an inflection point whose the 
first angle difference is an even multiple of tt and Its 
end points are at the steps whose first angle differ¬ 
ences are odd multiples of tt. This gives each con¬ 
jugate region the characteristic “f-slot” shape on the 
belly of a stringed instrument, to be described in de¬ 
tail below. 

Figure m the plot of f/rr ~ 300 million steps for 
s = 1/2 + 10®i, illustrates these details. For this 
large value of t, the angle difference regularities of 
alternating odd and even multiples of tt are clearly 
evident for steps and their conjugate regions. A red 
dot marks the center of symmetry of this step plot, 
which shows that the direction angles of each Initial 
step and its conjugate region add to a constant, whose 
perpendicular defines an axis of symmetry (dashed in 
the figure). In particular, since the first step is at the 
angle 0 with length 1, the angle of the final conjugate 
region is this constant angle sum, denoted here by 
20(s). 

A central proposition of this work is that the geo¬ 
metric center of each scroll corresponds precisely to 
the conjugate point of an initial step. Because of 
step discretization, this differs from the midpoint of 
a step by a longitudinal displacement A which de¬ 
pends on a and a transverse displacement irft) from 
the step center. A first-order derivation of these dis¬ 
placements is given below. 

The density of steps near the symmetry center of 
this figure precludes visualization of Its details; these 
are more clearly displayed in Figure [21 which reveals 
other Interesting features such as the central pendant 

^Although derived from a different approach, this quantity 
plays a significant role in Riemann’s analysis. It can also be 
defined by the condition S^8 = 2it. 
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with its self-symmetric center and fractal-like repe¬ 
titions where angle changes reproduce smaller repre¬ 
sentations of the pendant. There are many conjugate 
figures having rational ratios of steps in their respec¬ 
tive counterparts. These details are described and 
illustrated in [7]. 

5 Extent of Conjugate Regions 

The extent of the region conjugate to step n < Up is 
given by 


N„-i 

Nr, = 

Nr, 

The inflection point of the region conjugate to step 
n is at step N° = [np/(n—1/2)], where the step angle 
is 9^0 = —tlog{N°), the first angle difference mod, 2 tt 
is 0, and the second angle difference is 27r(np/A^°)^. 
An approximate result for the sum is given by Ja¬ 
cobi’s 0-function G{u). This is appropriate for large 
t and n sufficiently less than Up that there are many 
steps in the sum. The approximately linear variation 
of step length across the conjugate region can be ne¬ 
glected, since it sums to zero, and the second order 
variation is Similarly, the third angle dif¬ 

ference is 0(N^)~^. To this order of approximation, 
step angles can be assumed to vary as the square of 
Sn = n — N°. This leads to the classic Cornu or 
Euler spiral, where the angle of the tangent to the 
curve varies quadratically with arc length. The lim¬ 
its of the sum can be taken as ±oo because the spirals 
near the end scrolls oscillate around a central value, 
giving 



'n- 1/2 


- 1/2]. (10) 


N, 


(n-l) 

E ^ 

Nn 




'E 


iTrSn^ —^ 


( 11 ) 


The approximatioiH G{u) = ® ~ 1/u 

gives 


N„- 


? = — - i 

Pfl-s 


i29^Q —i7r/4+.. 


Nrt 


^Valid to five places for |w| = 


( 12 ) 


where is the angle of the step at the symmetry 
center. Algebraically, the angle 7r/4 is the term \fi 
in the definition of u; from Figure [T] it is seen to be 
the angle between the central step about which the 
integral was evaluated and the summed extent of the 
Euler spiral. A further angle term, denoted by 5t, is 
due to the difference between t and the exact value for 
which 5^9no = ttN'^/ n^. Because the reference value 
of St is an integer multiple of 27r, this correction is de¬ 
noted in [3] simply as t. Higher order angle terms due 
to the approximations of this approach are bounded 
by the first omitted term from Riemann’s analysis, 
l/48t. 

The ratio of Ewf n ® to ® is central to this 
work, as it quantifies the overall symmetry of the 
step plot. It is denoted here by 


Q(S) = „l-2<Tg2*e„^-^/4+i+... ^ „l-2<Tg2*e (^3) 

which, except for the choice of sign, is identi¬ 
cal to Riemann’s result, as seen in the functional 
equation[31 Sec. 6.5]. It can be shown that ({s) 
lies on a line through the origin at the angle 0 = 
9np — -l- t/2 - 1 - .... 

The factor (5(s) can be related to the results of 
analytic number theory. Consider the comparison of 
the integral around the nth imaginary axis pole pair 
of the generating function to the contour integral of 
the single step n\ 

{2TTi){—2'iTiny~^ O r(s)n~^ (14) 

The integral of the nth pole pair on the imaginary 
axis corresponds to the sum of the steps conjugate 
to step n. Riemann used the T-function identities to 
express this ratio in a form emphasizing the symme¬ 
try under s —>■ 1 — s, as it is usually presented in the 
literature. 


C(s) = r(l — s)(27r)® ^2sm(s7r/2)^(l — s) cr > 1. 

(15) 

However, the ratio can also be expressed us¬ 
ing Stirling’s approximation r(s) = (27r)^/^(s — 
l)s-i/2gi-s ^ giving 7ji-2sgi(7i./4-i-t-i-...) (;je]-ived 
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Figure 1: Plot of summed steps for C(l/2 + 10®i). 



Figure 2: Steps near the symmetry center for s = 1/2 + 10®i 
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above. This form of Q{s) emphasizes the symme¬ 
try of initial steps and conjugate regions, specifically 
displaying the dependence of |(5| = on t. 

The functional equation of ^ is an immediate con¬ 
sequence of Equation [T^ except that this only de¬ 
scribes “half” of the relation, the steps and 
for n < Up. Equation [T5] results from the sum of 
Equation 1121 and its transform under s —> 1 ~ s, each 
summed to rip, and the relation Q{s)Q{l — s) = 1. 


6 First-Order Determination of 
the Conjugate Points 


Figure [T] strongly suggests that the point conjugate 
to the end of step n is the center of the scroll having 
56 = — (2n -I- l)7r/2, or t = {2n + l)7r. For this exact 
value of t, second order evaluation of 50]y-i/2±i/2 
gives the symmetric trio of steps iV — 1, TV, TV -|- 1 
with equally spaced angles tt 7r/2TV. When t equals 
this odd multiple of tt, the scroll center lies on step 
TV, giving a transverse displacement tat = 0. The 
midpoints of steps TV, TV -|-1 are separated by one half 
the step length TV”'* times their angular separation 
7r/2TV. Defining 5t = t — 7r[t/7r], each increase of St 
by TT advances the scroll center to subsequent steps 
TV -I- 1, TV -I- 2, etc. Although the angle changes of the 
steps are large, equal to 6Uog{N), the variation of 
the angle difference 56 n with 5t is small, yielding the 
linear relation tn = T(5t/4TV®+^. 

For cr ^ 0, the variation of step length with n 
leads to a longitudinal displacement Ajv between the 
midpoint of a step and the geometric center of the 
scroll. Defining that center as a fraction a of its 
step length, Aat is given by the condition that the 
centers of adjacent steps having angle changes of tt 
coincide. Assuming that a varies weakly with n, 
(1 — a)TV“® = a(TV -|- 1)“® gives a = 1/2 -|- ajS^N. 
The effect of this longitudinal displacement on Tjq is 
of higher order. 

Taken together, and including the overall step an¬ 
gle —tlog(n), the scroll center conjugate to step n = 
T/tt is given to first order by 


N 


E 


m 


1 


1 

2N^ 


a + iSt 
4TV^ 


(16) 


The inclusion of higher order corrections for com¬ 
parison with Riemann’s analytic results is left to the 
reader. 


7 A Point Conjugate to Itself 

The complete step plot has bilateral symmetry, but 
in general the point Up does not lie on the symmetry 
axis. Figured] clearly shows the “pendant” near step 
Up and its polygonal shape which approximates a cir¬ 
cle for Up^ 1. The center of this pendant is a point 
on the axis, denoted here by P{s). It is laterally dis¬ 
placed from X)!*” by the distance L, as shown in 

Fig. El 



Figure 3: Location of the Symmetry Center 

Because 5^6 = 27r at step Up, the steps Up and 
Up + l are tangent to this central feature. Geometrical 
construction gives 


Tip 

P{s)=^n-^ + L (17) 

1 

where 
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(18) 


7V-" 


g-i(tiog(rtp)+27rp) 

2np cos(27rp) 


and where p is the fractional part of •\/t/27r, given 
by {np + p)'^ = t/2TT. Note that L 0, as seen both 
by the construction and Equation [T51 Figure [3] also 
shows that L is closely parallel to the axis of symme¬ 
try, as it is perpendicular to the average of angles 
and 9n+i- 


N 

c(s) = 

1 


s - 1 ^ (2fc)!iV2fc ■ 


(19) 

where the Bernoulli numbers Bn are given by the 
generating function 



1~) 

1 


nl 


( 20 ) 


8 The Point Conjugate to the 
Origin 

The reciprocity relation between original and con¬ 
jugate steps implies a final scroll near step [t/7r]. 
Consider a pair of adjacent steps n,n + 1 beyond 
the inflection point of the final conjugate region at 
</27r, where SO decreases uniformly to zero with no 
further discontinuities of 2tt due to the mod function. 
The perpendicular bisectors of these adjacent steps 
intersect at a distance R = n~^jSQ ~ jt, 
defining a radius of curvature that decreases with 
increasing n for a > 1. In this case, the final scroll 
collapses to the value of When cr < 1, i? increases 
without bound. 

The series divergence is entirely due to the 
outward spiral of steps beyond the final scroll at 
O' = cr < 1. 


The sum over powers of {s/N) converges slowly 
for |s/iV| >= TT, but it can be represented by the 
Bernoulli number generating function. In this way, 
one obtains the first order correction of the final step; 

jsj-s g ^s/2N^^-s/2N 

_g-s/2N ) ) 

The ratio s/(s — I) approaches I for large s, and 
manipulation of the definitions of trigonometric and 
hyperbolic functions lead to 

tanh(cr/27V) — i ■ cot{t/2N) 

2 1 — i • tanh(cr/2iV) cot(t/27V) 

Expansion of the real term reproduces the longi¬ 
tudinal correction term cr/4A^'*+^ for the final scroll. 
The term cot{t/2N) implements the mod function, 
and its expansion about the zero of the cot gives 
resulting in agreement with the geomet¬ 
ric first order evaluation of the point O' conjugate to 
the origin. 


Edwards notes that Riemann did not use the term 
“analytic continuation”, but rather discussed “a so¬ 
lution that remains valid” when the series representa¬ 
tion does not converge. By that definition, the point 
O' conjugate to the origin is the analytic continua¬ 
tion of C(s) for cr < 1. This conclusion appears to 
differ from all methods reported previously, but is in 
agreement with them. As one example, consider the 
Euler-Maclaurin summation method, in which the ex¬ 
plicit sum of terms beyond n = [t/27r] is represented 
by an integration which extends the region of con¬ 
vergence to cr > 0. This yields, to lowest order [3l 
Sec. 6.4] 


N=[t/pi\ 


1 


N ^ a -\-iSt 
^ 4iV«+i 


(23) 


9 First-Order Riemann-Siegel 
Correction 


Although Riemann did not include it in his paper, his 
handwritten notes (recovered by Siegel in a feat of 
exceptional mathematical ability) show that he used 
symmetry to facilitate calculations, including a hi¬ 
erarchy of equations giving successively higher order 
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results. His first order result can be shown to corre¬ 
spond to 2L{s), the lowest order of Equation [T8l 

In Riemann’s development, its conjugate 

Q{s) and a remainder term R bridging the 

endpoints of the two sums are projected to the 0-axis 
on which C(s) resides, and termed Z(s). , 

Tip Tip 

1 1 

Since Z S 5R, locating its zeros is a matter of de¬ 
tecting sign changes. Riemann determined the first 
order “mismatch” between s-nd its conju¬ 

gate using the “steepest descent” or “saddle-point” 
method[3] Sec. 7.2]. The integrand of the mismatch 
is expanded about its peak. One segment of a closed 
contour crosses the imaginary axis at the Upth pole 
at the angle 7r/4, leading to an integral of gaussian 
form with its peak on the imaginary axis. All other 
portions of this contour give negligible contributions. 
His first order expression, including his approxima¬ 
tion for the remainder R and assuming tr = 1/2, is 

n-p 

Z = n“^/^2cos(0(t) -I- tlog{n))+ 

1 

( 25 ) 

Up cos(27rp)) 

From symmetry, L is related geometrically to Rie¬ 
mann’s remainder R by the expression R = L(s) + 
Q{s)L{l — s). In the approximation that L lies ap¬ 
proximately along the symmetry axis this can be sim¬ 
plified to 

R = 2Lcos(6»l - 0) (26) 

where 

0 = —tlogijip) + TTUp + 7r(— — 2p^) 

TT 1 

6*l = -tlog{np) - - - -66 

(27) 


Using Equation [TH| and noting that 0 involves p 
through tlog{^yjt/2TT)) = tlogijip) -|- tlog{l + p/up), 
discarding terms of 0{l/np), and using Trnp{mod2Tr) 
gives 0 for Up even or tt for Up odd gives 

0L — 0 = 27r(p^ — p — 1/16-|-(0 or tt)) (28) 

Combined with the equation for L, this shows that 
Riemann’s first order correction is given by the pro¬ 
jection of 2L to the 0-axis. A study of Riemann’s 
analysis in [31 Sec. 7.4] displays both his mathemat¬ 
ical proficiency and the simplicity of the geometric 
approach. 

The zeros in both the numerator and denominator 
of i? at p = 1/4, 3/4 can be understood geometrically. 
At these values of p, 6dnp is a multiple of 27r and 
steps rip and Up+i are each exactly perpendicular to 
the symmetry axis. The error in the pendant center 
location grows large at these values of p, but does so 
along the axis. It is therefore perpendicular to 6 l. In 
the projection to the 0-axis, the zeros of cos(0 — 6^) 
cancel this effect. 

10 The “Double Pendulum” 
Construction 

The Riemann-Siegel calculation is represented geo¬ 
metrically by the sum of two line segments, OP — 
P{s) and O'P = Q{s)P{l — s). This “double pendu¬ 
lum”, related to a classical geometrical figure called 
the limagon, is illustrated in Figured] for s = 1/2 -|- 
2220000.15i, chosen to reduce conjugate region over¬ 
laps and display the features clearly. The green vec¬ 
tor connects the origin to the pendant center P{s). 
This, plus its conjugate, in black, ends at C(s)- The 
bisector of these two lines is the symmetry axis, show¬ 
ing geometrically that ({s) lies on the line through 
the origin perpendicular to the symmetry axis. 

Figure [5] illustrates this construction for a range 
of values of t having a = 1/2. Here, the dashed 
spiral shows P{s) and the solid spiral gives O' = ({s) 
at successive values of t, with the curves connected 
at regularly spaced values of t to indicate the line 
segments of the limagon construction. 
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Figure 4: Riemann-Siegel illustration for s = 1/2 + 
2220000.15i. 



Figure 5: Limagon construction for Gram points 
1000 ^ 1004. 



Figure 6: Distribution of fractional positions of the 
zeros between adjacent Gram points, for the first 
100,000 zeros. 


The number of zeros for t < T has been a subject of 
much research. Integration around the boundary of 
the critical strip shows that there are approximately 
as many zeros as the number of times that 0(T) has 
passed through a multiple of n. This leads to the 
definition of Gram points hy the relation ©(Itv) = 
IVtt [3J Sec. 6.5]. When 0 is an odd multiple of 7r/2 
the real part of ( is zero, countering its “bias” to 
positive values because the first step always gives a 
real contribution of +1. In this case, nearby values of 
t are more likely to lead to complex zeros. Although 
the zeros are spaced irregularly because the angle of 
P{s) does not change monotonically with t as does 0, 
on average there does appear to be one zero per Gram 
point and they are found more frequently between 
successive Gram points. The tendency for zeros to 
“group” between Gram points is shown in Figure [6l 
a histogram of the displacement of the first 100,000 
zeros [8] from the midpoint of adjacent Gram points, 
scaled by half the distance between Gram points. 

When a = 1/2, the condition of opposing angles is 
a one-parameter “search” as t increases, just as the 
sign changes of Z signal the crossing of a complex 
zero for C{t). Thus, the symmetry shows why zeros 
occur frequently for this value of a. This has led to 
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Figure 7: 
Q{s)Pil - 
t < 129. 


Magnitudes of P{s) (plain mesh) and 
s) (color) on the critical strip for 124 < 


Figure 8: Euler-Maclaurin, 2000 < t < 2010, blue : 
a = .500, red : a = .505 


the Riemann hypothesis [H Sec. 13.9] that all zeros 
are on this line in the complex plane. Stated in this 
form, a proof of the hypothesis requires showing that 
zeros are impossible when cr ^ 1/2. As illustrated in 
Figure [3 P(s) and Q{s)P{l — s) for 120 < t < 129 
intersect transversely in this portion of the critical 
strip. This intersection is identical for a = 1/2, where 
the transformation s —>■ 1 — s corresponds to complex 
conjugation and the magnitude of Q{s) is unity. 


11 Zeros and the Critical Line 

Riemann’s On the Number of Prime Numbers Less 
Than a Given Quantity was chosen for inclusion in 
God Greated the Integers^, Hawking’s list of land¬ 
mark mathematical works. Riemann found that the 
first few zeros are on the critical line. However, this 
result was not obvious from the calculational algo¬ 
rithm; this peculiar generality was later stated as the 
general hypothesis bearing his name. Recognition of 
the symmetry sheds light on this, as can be illustrated 
by three specific algorithms for f. Since each is a 
solution of the same problem, presented from differ¬ 
ent viewpoints, they give the same numerical results 
within their respective error limits. 

The first algorithm evaluates C, for a range of t and 


arbitrary values of a using the Euler-Maclaurin ap¬ 
proximation; 


N—[t/-n\ 

C(s) ^ W 
1 


1 


a + iSt 
4_/V®+i 


(29) 


Figure |5| shows versus JJC for 2000 < t < 2010, 
with a = .500 in blue and a = .505 in red. En¬ 
largement of the region around the origin in Figure 
El shows that every “loop” with cr = 1/2 invariably 
yields a complex zero, whereas a slight variation of a 
gives “loops” which avoid the origin. 

How can this extremely simple algorithm give a 
zero for every “loop”, but only when tr = 1/2? When 
the author first encountered this result, having then 
no appreciation of the symmetry, it seemed wildly 
unusual, being nothing more than the sum of Dirich- 
let terms with the first order Euler-Maclaurin 
correction. What “unseen hand” could possibly force 
this particular result? With symmetry, one realizes 
that although P{s) and O' are not explicitly rep¬ 
resented in in this prescription, their effects are in¬ 
cluded correctly in the result. Note that the Euler- 
Maclaurin algorithm is insensitive to the upper limit 
N > [t/27r]; terminating the sum at iV -|- 1 and sub¬ 
tracting {N -I-1)“'’/2 gives the center of the next step. 
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Figure 9: Euler-Maclaurin enlarged, 2000 < t < 
2010, blue : a = .500, red : a = .505 

and the change of St by tt yields the same value for 
the spiral center as did the upper limit N. 

Use of Riemann’s formulation in Equation [25] gives 
the same result far more efficiently, and in compact 
notation, summing only 17 terms instead of more 
than 600 for the Euler-Maclaurin algorithm with 
same range of values of t. 

tip 

Z = n“^/^2cos(0(t) -I- tlog{n))+ 

1 

po) 

Up cos(27rp)) 

However, this form is limited to a = 1/2 due to 
the compact expression of the projections of 
and to the 0-axis using the cosine function. 

The remainder term is presented as a single entity 
bridging the two discrete sums, masking the identity 
of the symmetry center. 

Implementation of the symmetric form 

C(s)=P(s) + Q(s)P(l-s) (31) 

using Equations ITTI and [TSl demonstrates the symme¬ 
try explicitly and is valid for arbitrary a. It gives the 



Figure 10: C Real (-) and imaginary (- .) components 
vs. t for Riemann’s form (blue) and the symmetric 
form (red). Red diamonds are the actual zeros. 

results and numerical efficiency of Riemann’s form, 
as shown in Figure 1101 Seen from this perspective, 
however, there is no peculiarity in the preference for 
zeros when a = 1/2. 

The functional equation demands that zeros off of 
the critical line must appear in pairs ats=l/2±e-|- 
it, and exceptional conditions are required for this 
to occur. The symmetric Riemann-Siegel equation 
shows that either: 

The magnitudes of P{s) and Q{s)P{l — s) 
must be equal for zeros at a = 1/2 ± e 
and the condition of opposing angles must 
be met, 

or: 

Since Q(s) ^ 0 on the critical strip, P{l/2 + 
e + it) and P(l/2 — e + it) must be simulta¬ 
neously zero at these arguments. 

Symmetry, with the example calculation shown in 
Figure|3 exposes the difficulty of attaining these con¬ 
ditions; it may, however, be the case that they cannot 
be excluded “in the fullness of t”. If so, this would ex¬ 
plain why the Riemann Hypothesis has resisted any 
proof for over 150 years [5]. 
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A detailed, internal symmetry exists between individual terms 
where n G P is less than a particular value rip, and sums over 
conjugate regions consisting of adjoining steps n greater than Up. 
The argument s — cr it is generally assumed here to lie in the crit¬ 
ical strip cr E (0,1), with t >> 1. Terms n“®, referred to as steps, 
are described by their lengths n~‘^ and angles On — —tlog{n). The 
boundaries of the conjugate regions are where first angle differences 
69n — —tlog{n + 1/n) equal odd multiples of tt. Two significant points 
in the complex plane are defined by this symmetry: 0'{s), conjugate 
to the origin O, and which equals C(s) not only for cr > 1 but also for 
for cr E (0,1); and P{s), conjugate to itself, which gives Riemann’s 
correction to the discrete sum in the Riemann-Siegel equation. The 
distances from P to O and P to O' are equal only for cr — 1/2. Su¬ 
perposition of O and O' results under the single-parameter condition 
that /LOP and ZPO' are opposed, giving frequent, although irregu¬ 
larly spaced, zeros of Ci^)- Analysis of this symmetry allows an alter¬ 
nate understanding of many of the results of number theory relating 
to C('S): including its functional equation, analytic continuation, the 
Riemann-Siegel equation, and its zeros. Discussion of three explicit 
computational algorithms illustrates that the apparent peculiarity of 
the occurrence of zeros when cr = 1/2 is removed by direct recogni¬ 
tion of the symmetry. The difficult conditions required for a zero 
when cr 1/2 are presented. 
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